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Abstract 

We investigate the validity of the Bogoliubov c-number approximation in the case of 
■ interacting Bose-gas in a homogeneous random media. To take into account the possible 

occurence of type III generalized Bose-Einstein condensation (i.e. the occurrence of 
condensation in an infinitesimal band of low kinetic energy modes without macroscopic 
occupation of any of them) we generalize the c-number substitution procedure to this 
band of modes with low momentum. We show that, as in the case of the one-mode 
condensation for translation-invariant interacting systems, this procedure has no effect 
on the exact value of the pressure in the thermodynamic limit, assuming that the c- 
numbers are chosen according to a suitable variational principle. We then discuss the 
relation between these c-numbers and the (total) density of the condensate. 

d ' 

Keywords: Generalized Bose-Einstein Condensation, Random Potentials, Bogoliubov 
c-number Approximation, Bogoliubov Quasi- Averages, Berezin-Lieb Inequalities. 

X 



in 

(N 



PACS: 05.30. Jp, 03.75.Fi, 67.40.-w 
AMS: 82B10, 82B23, 81V70 



1 PhD student at UCD (Dublin, Ireland) and Universite de la Mediterranee (Aix-Marseille II, France), 
e-mail: Thomas. Jaeck@ucdconnect.ie, phone: +353 1 7162571 

2 e-mail: Valentin.Zagrebnov@cpt.univ-mrs.fr, phone: +33 491 26 95 04 



1 



1 Motivation 



In 1947 Bogoliubov pQ proposed an ansatz that for large Boson systems the particle cre- 
ation and annihilation operators: a^, clq, corresponding to the zero momentum k — 0, can 
be replaced by complex numbers. This procedure is known as the Bogoliubov c-number 
approximation. It is based on the idea that these creation and annihilation operators, when 
divided by the square root of the volume, V, of the region A containing the system, can be 
expressed as space averages: 



# 



V 



i j dx a*(x), 



where aft(x) are the usual local creation and annihilation operators. Therefore for translation- 
invariant ergodic states these operators should converge to multiples of the identity: 



a* 



V 

in some weak sense, see e.g. [3]. These ideas were exploited to construct a various truncations 
of the full interacting boson Hamiltonian. We refer the reader to, e.g., [2] for a review of 
these models and to [22], [23] for some recent applications. 

The most spectacular result derived from this ansatz was its application to a homogeneous 
model of a weakly interacting Bose gas [1], [13], which provided explicitly a spectrum of col- 
lective excitations satisfying the Landau criteria of superfluidity. Note that this microscopic 
theory of superfluidity is also based on two other Bogoluibov's ansatse: the occurrence in 
the weakly interacting boson system of condensation in the k = mode and the truncation 
of the full Hamiltonian, keeping only the "dominant" terms, that is those that involves at 
least two particles from the condensate. 

Recently, see [IT], [T8], the Bogoliubov approximation has been used to study interacting 
bosons systems in homogeneous external random potentials, where the notion of the ground- 
state as well as the existence of the condensation are quite subtle. The aim of our paper is 
to study the validity of the Bogoliubov scheme for this kind of models. 

The first rigorous result concerning the Bogoliubov c-number approximation was due to 
Ginibre [I] . For a homogeneous boson gas with a two-body superstable interaction he proved 
that the Bogoliubov ansatz, supplemented by a self-consistency condensate equation which 
maximizes the approximated pressure with respect to the c-number a, gives the right pres- 
sure in the thermodynamic limit. A transparent and elegant proof of this and other related 
results has been recently given by Lieb et al [5] , using, in particular, the Berezin-Lieb and the 
Bogoliubov convexity inequalities. The paper [5] also investigates a delicate point: namely, 
whether the value of the variational parameter a max maximizing the approximating pres- 
sure coincides with the total condensate density. There it was shown that the maximizer 
« max corresponds to the zero-mode condensate density if the gauge symmetry breaking term 
(quasi- average sources) of the form: \/V(r]aQ + fja ) is added to the full Hamiltonian. The 
idea of breaking the gauge-symmetry of the quantum Gibbs state for k = is due to Bo- 
goliubov, [2J . This forces the totality of the condensate to be concentrated in the zero-mode 
(ground state). It must be switched off (rj = \r)\ e iargri , \rj\ — >• with a fixed gauge <fi := argr/) 
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after the thermodynamic limit to produce a limiting Gibbs state. The expectation defined 
by this state is called the Bogoliubov quasi- average with respect to this source. The quasi- 
averages of the operators af /y/V coincide with |a max | e ±l< ^, i.e., |a max | 2 is equal to the total 
condensate density [5]. It has been argued in [2], j5] that this quasi-average is the only 
physically reliable quantity to describe Bose Einstein condensation. 

We emphasize this point, because it has been known since [6] that the Bose-Einstein 
condensation in the gauge-invariant systems does not necessarily imply a macroscopic occu- 
pation of the ground state only. Indeed, the condensate can be spread over many (possibly 
infinitely many) quantum states, and in some cases, none of these states are macroscopically 
occupied. In all cases however, the total amount of condensate in an arbitrary small band of 
energy in the vicinity of the ground-state is the same, a phenomenon known as generalized 
condensation in the terminology of van den Berg-Lewis-Pule [6]. 

After Ginibre [I] it was tempting to conjecture that the limiting value of the solution 
to the finite volume condensate equation yields the correct condensate density under all 
circumstances. By means of a counter example it was shown in [7J that this is not so. 
Although the Bogoliubov c-number approximation still gives the right pressure in these 
systems, it has been shown for the mean- field boson gas [TJ, that the solution a max of the 
condensate equation does not provide the ground-state condensate density po, but the total 
amount of the generalized condensate, i.e., po < |cw x | 2 - This was in a striking contrast with 
a general conviction that the gauge-invariant, translation invariant boson systems always 
manifest the total amount condensation in the ground state, and hence that the Bogoliubov 
c-number approximation in the zero- mode coincides with this amount. 

We return to this point at the end of the present paper where we discuss our main result 
(Theorem 13.11) in relation to the results of [I], [TJ, [2], and [5j. 

The aim of our paper is to prove that the Bogoliubov c-number approximation to inter- 
acting Bose-gas can be extended to the case of homogeneous random external potentials 
(random media). Note that the arguments in [5J are valid for inhomogeneous systems and 
allows to treat condensation in many modes, as long as its number is much smaller than 
the volume V. In the case of random media this problem is more complicated. First, we 
have to investigate a random inhomogeneous system, albeit with non-random properties 
(self- averaging) in the thermodynamic limit [8]. Secondly, since the randomness seems to 
force the generalised condensation to be of type III, as we disucussed in [15] . the number of 
modes occupied by the condensation in the random potentials is a priori of the order of V. 

Recall that for non-interacting (perfect) bosons systems embedded into a bona fide ran- 
dom potential, the generalized condensation occurs even in low- dimensional systems (see 
0, QUI E] and [12], [15] for a review) even though it does not occur for the corresponding 
translation-invariant systems. This is caused by the fact that the one-particle density of 
states has Lifshitz tails, that is an extremely low density of quantum states near the bottom 
of the spectrum, a well-known feature of random systems widely believed to be associated 
with the existence of localised eigenstates. Kac and Luttinger [9j [10] conjectured that in 
a homogeneous random potential condensation occurs in the state with the lowest energy 
(ground state). Indeed, one can check this conjecture [13] for the particular case of the 
Luttinger-Sy model [TT]. Recently in [T5J, it was shown that whenever condensation occurs 
in the random perfect Bose-gas, then there is a generalized condensation in the kinetic-energy 
momentum states, and the both densities have the same value. This result can be partially 
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extended to some simple models of interacting Bose gases (mean- field models). It was proved 
in [IB] that under a fairly weak assumption about localization this condensation is of the 
type III, i.e. there is no macroscopic accumulation of particles in any single momentum 
state. We conjecture that this holds also for general interacting systems in the presence of 
homogeneous random potentials. 

For the above reasons, we would like to follow the general philosophy of the Bogoliubov 
c-number substitution ansatz and extend it to a generalized Bogoliubov approximation, in 
order to cover the possible case of type III generalized condensation. By this, we mean a 
replacement of all creation/annihilation operators corresponding to momentum states with 
kinetic energy Sk in the energy band < Sk < 8 by complex numbers {VV &k}{k-e k <8}- 
First we show that this extension of the Bogoliubov approximation applied to interacting 
Bose-gas in homogeneous random potentials is valid as far as the pressure is concerned 
(Section [3]). In this case as in the case one- mode condensation the corresponding trial 
pressure is maximized with respect to these complex numbers and then one lets the parameter 
5—7-0 after the thermodynamic limit. Note that for each realization of the random potential, 
the system is not translation invariant and the minimizer has a random value before the 
thermodynamic limit. For this reason, the proof use the stationarity and ergodicity of the 
random potential with respect to the space translations. 

Finally, we discuss the variational problem established for the pressure (Section |4]). In 
particular, we highlight the fact that the link between the c-numbers that maximize the trial 
pressure and the structure of the condensate is highly non-trivial. By the mean of a simple 
example, we show that the Bogoliubov quasi-average technique of adding external sources [2] 
is not "satisfactory" in this case, since we suspect that the generalised condensate should be 
of type III, while this procedure is able to drastically alter the fine structure of the generalised 
condensate. This brings us back to the discussion of the gauge symmetry breaking and the 
physical reliability of the Bogoliubov quasi-averages in the case of condensation [2], 0. 



2 Model and definitions 

Let {A; := (— 1/2,1 /2) d }i^i be a sequence of hypercubes of side / in M. d , al ^ 1, centered at 
the origin of coordinates with volumes V% = l d . We consider a system of identical bosons, 
of mass m, contained in A/. For simplicity, we use a system of units such that h = m = 1. 
First we define the self-adjoint one-particle kinetic-energy operator of our system by: 

^:=-IA p , (2.1) 

acting in the Hilbert space M{ := L 2 (Ai). The subscript P stands for the periodic boundary 
condition. We denote by {ip l k , e l k }keA* the set of normalized eigenfunctions and eigenvalues 
corresponding to operator h® 

H(x) = 4 = \k\ (2.2) 

and A ; * is the usual dual space to A;, that is A^* := {k G M. d : k 2 = (27m) 2 /l 2 ,n G I* d }. 
Finally, we denote by uf the finite- volume integrated density of states (IDS), that is, 

v?(E) ■= il{^A?: £ U4 (2.3) 
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and we let u°(E) := lini/^oo uf(E). Note that the limiting IDS, u°(E), has support on [0, oo) 
and that it is known explicitly: u°(E) = CdE d l 2 (the Weyl formula). 

Definition 2.1 We define an external random potential v^'(-) : f2 x M. d — >■ K, 14 as 
a measurable random field on a probability space (f2, J 7 , P), satisfying the following conditions: 

(i) G f2, non-negative; 

(ii) p := P{w : ^(0) = 0} < 1. 

As usual, we assume also (see e.g. Appendix B in [15]) that this field is: 

(hi) homogeneous (stationary) and ergodic with respect to the group {r x } xe ^d of probability 

preserving translations on (Q, J 7 , P) ; 

(iv) (p-mixing for S ^-measurable functions, where £a is the a -algebra generated by the field 
{v u (x)} xeA for Ac R d . 

Then the corresponding self-adjoint random Schrodinger operator acting in J$f := L 2 (M. d ) is 
a perturbation of the kinetic-energy operator: 

h u := -\ A + v", (2.4) 

defined as a sum in the quadratic-forms sense. The restriction to the box A/, is specified by 
the periodic boundary conditions and for regular potentials one gets the self-adjoint operator: 

hf:=H A + <) p = (2.5) 

acting in Jfi, where vf is the restriction of v w to A/. We denote by {4>f'\ Ef'}^ the set of 
normalized eigenfunctions and the corresponding eigenvalues of the random operator hf. We 
order the eigenvalues (counting the multiplicity) in such a way that Ei sC E% ^ K 
Note that the non-negativity of the random potential implies that Ef' 1 > 0. So, for con- 
venience we assume also that in the thermodynamic limit the lowest edge of this random 
one-particle spectrum o~(hf) satisfies the fifth condition: 

(v) lim^oo Ef' 1 = , almost surely (a.s.) with respect to the probability P. 

Remark 2.1 Note that (v) is in fact an implicit condition on the random potential saying 
that a.s., one can find a sequence of regions (gaps) with v w (x) = 0, with volumes tending to 
infinity in the van Hove sense. 

Now, we turn to the many-body problem. Let &\ := be the symmetric Fock space 

constructed over 3%\. Then H® := dT(hf) denotes the second quantization of the one- 
particle Schrodinger operator hf in For simplicity, we omit the explicit mention of the 
randomness of the Hamiltonians and all related quantities, unless this is necessary for the 
sake of clarity. Since for any oj G VL the one-particle eigenstates {<fri := 0"' }i>i of hf form a 
basis in j£f, the operator Hf acting in &\ can be expressed as: 

Hf = = E^'^W • (2-6) 

Here a*(0j), a(4>i) are the creation and annihilation operators, satisfying the boson Canonical 
Commutation Relations, and iVj(0j) is the particle-number operator in the state <pf' 1 . Note 
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that, since [h®, hf] 7^ 0, the Hamiltonian (12. 6p cannot be expressed as a function of the 
operators Ni(ip l k ) in the kinetic-energy eigenstates (12. 2p . 

(vi) Below we assume that the particles interact through a suitable non-negative two-body 
translation-invariant potential $(x, y) := u(\x — y\). More precisely, we assume that the 
function u has a continuous, bounded Fourier transformation u(q), and there is 7 < 00 such 
that \u(q)\ < 7 for all q G A z * and for all I. For example, one can choose u e L 1 (M rf ). 

Then the second quantization of interaction: C/j := dr($) has a simple form in the 
translation-invariant basis {ifjk '■= V^lfceA* °f the kinetic-energy operator h®: 

Ul = ~jy Yl ^ a *(i>k+q) a *(i J k>- q )a( y i} k >)a(^ k ) . (2.7) 

1 q,k,k'eA* 

The full Hamiltonian with the chemical potential \x included has the form: 

Hiifi) := Hf-nty + Ui (2.8) 

Note that the creation and annihilation operators in the interaction term (12. 7p are in the 
momentum eigenstates {ipk}k£\*, although the perfect Bose-gas Hamiltonian (12. 6p is not 
diagonal when expressed in this basis. 

By (-)j we denote below the grand-canonical equilibrium state defined by the Hamiltonian 

(A)i(P,») := ^^Tr^exp(-/^(//)), 

and by pi({3,/i) its associated grand-canonical pressure 

p,0M := -LhiS^), (2.9) 
PVl 

where 

Ei{p,ii) := Tr^ ; exp(-/3^(/i)) 
is the corresponding grand-canonical partition function. 

It is known that the pressure of the corresponding non-random model (i.e. for v w (%) = 0) 
with a bona fide interaction exists and is independent of the boundaries condition for a large 
class of them, including the periodic case, see e.g. [19]. The proof of this statement consists 
essentially in showing the existence of the Dirichlet pressure using sub-additivity 

where A', A" are disjoints subsets of A, and r x denotes translation by x. The exact value of 
x is chosen according to the usual tempering condition required of the two-body interaction 
potential (vi). Then, using translation invariance of the non-random model, one obtains 

Pa CM > Pa'0M+Pa»(M , (2-10) 

and the boundeness of the pressure, which is provided by the super stability of the interaction 
{u > 0, (vi)), leads to the existence and finiteness of the limiting pressure for any fj,. Then, 
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one can show using functional integration techniques, see [20], that the others boundary 
conditions converge to the same limit. 

The last part of this prove can be carried through verbatim in the presence of the external 
random potential. However, because of the lack of translation invariance in the random case, 
the inequality (I2.10p for the Dirichlet pressure is modified as follows: 

Pa^^iu) ^ p% u (P,[i) +p°%,(P,fj) = p^{P,n) + p°f xUJ (f3,p) . (2.11) 

We used the stationarity of the random potential in the last identity. To prove the the 
existence of the thermodynamic limit one can use the Kingman sub-additive ergodic theorem, 
see |2T]: 

Proposition 2.1 Let r be measure preserving transformation of the probability space (fl, J 7 , P) 
and {g n }n>i be a sequence of functions g n G L l (Q, J 7 , P) satisfying the condition: 

g n + m {u) < g n {u) + g m {r n u) (2.12) 

Then one gets that 

a.s. — lim g n {uj)/n = g(u) , (2-13) 

n— ¥oo 

where the function g{oj) is r-invariant: g[r s bS) = g(u>). If in addtion, the functions g n are 
ergodic, it follows that the limit g{u) is a.s. non random. 



3 Generalized Bogoliubov c-numbers approximation 

3.1 Existence of the approximating pressure 

Following the Bogoliubov approximation philosophy, we want to replace all creation/annihilation 
operators in the momentum states ip k with kinetic energy less than some 5 > by c-numbers. 
To this end let Is C A ; * be the set of all replaceable modes, that is, 

I s ■= [k e A; : k 2 /2 <: 5}, 

and we denote rig := §{k : k G Is}- The number of quantum states rig is of the order V\, since 
by definition of the IDS ( 12. 3ft : ns = V% v®($). Let J^f 5 to be the subspace of M{ spanned 
by the set of ijr k with k G Is, and Ps be orthogonal projector onto this subspace. Hence, we 
have a natural decomposition of the total space <%{ and the corresponding representation for 
the associated symmetrised Fock space: 

M{ = ^®^, J^J£/®J^. (3.1) 

Then we proceed to the Bogoliubov substitution a k — > Ck and a* k — > Ck for all k G Is, which 
provides an approximating Hamiltonian that we denote by H\ ow ([jl, {c^}). The meaning of 
the superscript low will become clear in the next section. We postpone to Appendix A a 
description of the explicit form of this operator. The partition function and the corresponding 
pressure for this approximating Hamiltonian have the form: 

SfVfo}) := Tr^e-^° w ^», (3.2) 

P l i7(^{ck}) := hnZ l r(ii,{c k }) . (3.3) 
Vi 
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The principal result of the present paper is the following main theorem: 

Theorem 3.1 The c-numbers substitution for all operators in the energy-band 1$ does not 
affect the original pressure i2.9\) in the following sense: 




Note that the number of the substituted modes is of order V, since we let 5 I after the 
thermodynamic limit. 

3.2 Proof of the main Theorem 

Our method is a generalisation of the one invented in [5] . For the convenience of the reader, 
we postpone the proof of some technical lemmas to the next section. 

First we define a family of normalized coherent vectors in the Fock space with the 
vacuum state |0): 



labeled by the set of complex numbers {ck}kei 5 - With the help of vectors (13.51) . we define 
the lower symbol A low for any operator A in #j by the partial inner product: 



which is an operator in J^" 1 , see (13. ip . Similarly, using the tensor structure of (13.11) . we can 
define on ^ the upper symbol A up of the operator A in &\ through the integral represen- 
tation 



Here d 2 Cj := dRe(cj)dlm(cj)/n and \c)(c\ := (£) keIs \ck){ c k\ is the projector on the co- 
herent vectors (13.51) . Similar to the one- mode case one has the completeness property 
f £ns d 2 ci . . . d 2 c ns \c)(c\ = I. Note that, contrary to the lower symbols, the upper sym- 
bol does not necessarily exists, and it may not be unique. Although, this poses no problem 
in our case, since the existence (though not the unicity) of the upper symbols follows from 
the fact that the Hamiltonian (12.81) is polynomial in creation/annihilation operators. We 
postpone the explicit expressions of these upper symbols to Appendix IA1 

We then define two approximating Hamiltonians, that we denote by Hj ow (fi, {ck}) and 
i/" p (/i, {ca,}). They are obtained by replacing all creations and annihilations operators a k by 
their lower, respectively upper, symbols. We refer the reader to Appendix [A] for the details 
and the explicit expressions. 

Note that i^ low (/i, {cfc}) is obtained simply by replacing all operators {a k }kei s by the corre- 
sponding complex numbers {c k }kei 5 - Formally it corresponds to the Hamiltonian obtained 
by the standard Bogoliubov approximation, that is why it appears in Theorem 13.11 



(g) e -|cfc| 2 /2+c fe a* 
k€l s 



" |o> , 



(3.5) 



A low ({c k }) := (c\A\c) 
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In a way similar to ( 13. 2p . (1 3 . 3 j) . one can define by {ck}) the partition function for 

the Hamiltonian if z up (/i, {c^}), and by $ p s (//, {c^}) the corresponding pressure. 

Finally, we denote by (—)i ow an d (— ) up the grand-canonical equilibrium states related to 
the following (integrated) partition functions: 



'low 



/ d 2 Cl . . . d 2 c n ,Tr^xe-^ low ^-^», (3.6) 



EJ*0*) := / d^-.d^W^W^, (3.7) 

and we denote the associated pressures by p l °™ (fi) , p^(fi) . 

Now one can mimic the arguments of [3] and extend them to the multi-mode projections 
\c)(c\ case to produce the Bogoliubov-Peierls and the Berezin-Lieb inequalities for (13.61) and 
(13. 7p . These inequalities yield respectively lower and upper estimates for the grand partition 
function: 

d 2 Cl . . . d 2 c n , Tr^e-^° W W C *» < SjM < / d 2 Cl . . . d 2 c n , Ve"^^"- (3-8) 

C n s 1 Jc n s 1 

By a straightforward generalization of the arguments in [5] to the case of multi-mode coherent 
projection \c){c\, we obtain for any {ck}kei 6 the bound 

on the integrand in the left-hand side of (13.81) . This in particular implies 

max Tr^ e -^ ow W c *>) < H,(^) , (3.9) 

i.e. the estimate of the grand partition function from below. 

To find a similar bound for the right-hand side of (13. 8p from above, we note that if/ ow (/i, {c^}) 
and i?" p (/i, {c^}) are related by 

H? p ((i,{ck}) = ^ low (/i,{c fc }) + K(/i,{ Cfc }) , 

see Appendix 13 equation (IA.22I) for an explicit expression of k(/i, {c^}). If this is combined 
with the Bogoliubov convexity inequality: 



In / d 2 Cl . . . d 2 c ns Tr^e-^ P ^' {ck}) - In / d 2 Cl . . . d 2 c n ,Tr, 



e -m\ ov {^{c k }) 



fa, d 2 d . . . d 2 c n ,Tr^ ( - nfa { Cfc })e-^» p (*fa») 

then (13. 8 p and (13 . 10[) provide the upper bound: 

lnTr^e-^W < In / d 2 Cl . . . d 2 ; c na Tr ,±e~^ - (*(/*, {c k })) up . (3.11) 

JC n s 

Using the orthogonal projection P§ : M{ i— > J 1 ^ 6 , and in view of (1A.22[) . one can estimate the 
last term in (13. lip explicitly: 

-K(ji,{c k }) ^ Tr^hf-^Ps (3.12) 
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+ |(A + 2^(5) + 2^(25))^|c fc | 2 



kei s 

+ |(2^) + 2z/°(5))^4a fe 

fce/f 

< Tr((/# - fi)P t ) - 7^(5) (l - 4V^°(25) + |r/°(5) + V^°(2<^ 
+ 4 7 ^(25)(£(|^ 2 -l) + £a*a fc ) . 

Here I| := A ; * \ Then taking into account the upper symbol of the total number operator, 
we find that 

{<*}) + a( £(1*1 2 - 1) + £ a *°*) = ~ a ' ^ > ( 3 - 13 ) 

kei 5 keii 

which together with equation (13. Ill) provides the following estimate 

lnTr^e-^W ^ In / d 2 Cl . . . d 2 c ns Tr^e^ Hl rM^ (3.14) 

+ Tr((hf-p)P s )- 1 ^(6)(l-4V l uf(25) + ^uf(5)+V l uf(25)) 
+ 47^(2% In / d 2 Cl ...d 2 C „ 5 Tr^e-^ low ^» 

(3.15) 

To finish the proof, we need three lemmas. We postpone their proofs to the next section. 

Lemma 3.1 Suppose that the system of interacting bosons Ii2.8\) has a bounded limiting 
particle density p(fi) for any fixed ju£R: 

p(/i) := <9 M p(/i) := df, lim pi(/3,/S) < oo (3.16) 
see A2.9\) . Then one gets the estimates: 

lim sup lim sup -^—d^p 1 ™ (p) p(jj), (3.17) 

lim sup lim sup — d^p7 s {p) ^ p(ji) . (3.18) 

Next, we relate the integrated pressure p^ip) defined by (13. 6p to the maximum of the 
corresponding integrand. 

Lemma 3.2 For any a > 1, one has the estimate: 

In f d 2 Cl . . . dW c ' / ' fl,b ' M ( 3 - 19 ) 

Pvi Jc n s 

< ^lnmaxTr^e-^^" - -L ln(l - l/a) + ^Had,p\j {p)) 

PVI {c fe } < pV; /5 
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Notice that above statements are independent of the possible presence of random poten- 
tials. The next Lemma serves to control (ergodic) random external potentials. 



Lemma 3.3 Taking into account our assumptions on the random potentials, see Section® 
one gets the following inequality: 

lim sup ^Tr (h? - n)P s ^ u°(8) ((£-//) +E(v w (0)V) . 
Here E(-) denotes the expectations in the probability space (Q, J 7 , P). 

Returning back to the proof of Theorem 13.11 we get from (13.91) and (13.141) the estimates: 
max pgffo {c k }) < ftO*) < pgTO*) + ^r(Tr((^ - /*)P,) - 7^(5) (l - 4^(2*) + 
+ V^(2*))) + 4 7 ^(25)-U, In / d 2 Cl . . . dV&jrxe-^<***» 
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By Lemma I3.2[ this implies for any configuration cu the estimates: 



max p{«V, {c fc }) ^ Pi (jjl) ^ max p'J (jjl, {c k }) + K»(l, 8) , (3.20) 

{Cfe} {Cfe} 

where the random parameter K w (l, 8) is given by 

K u (l,8) : = ^(Trft/tf-^-T^^ 

+ 47^ (25)^ln / d 2 Cl ...d 2 c n5 Tr^e-^ low ^». 

Note that, by Lemmas 13.11 and 13 . 'S\ we can control this error term since for any configuration 
one gets: 



limliminf ^(1,8) = lim lim sup K"(l, 8) = . 

54,0 i— >-oo 54.0 ;_>.oo 



Hence, fl3~20]) yields 

limsup max {c^}) ^ p(fx) ^ limsup max p| ^(//, +liminf K u (l,8) 

l^oo { c k} ' l-^oo {cfc} ' 

limsup limsup max p'^"(//, { Cfe}) ^ p{ji) ^ liminf limsup max {c^}) , 

which proves the first equality in H3.4[) in Theorem 13.11 The second one can be proven in a 
similar way. 



3.3 Proofs of technical results 

Recall that the finite-volume IDS (12.31) : v®{8) > 0, converges in the limit Z — > 00 to the Weyl 
formula v°{8) = Cd8 d ^ 2 , uniformly in 8 on any finite interval < 8 < E. Below we use this 
uniformity in a systematic way. 
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Proof of Lemma 13.11 

Since 

Hj° w (fi, {c k }) +a(Y, N 2 + E a > k ) = H i™(» ~ a ' ^ > ( 3 - 21 ) 

kel s fcei| 

the estimate (13.121) yields: 

i/; ip (/i,{c fe }) ^ #^ + 4 7 ^(25),{c fe }) 

- Tr((h? - vl)P 5 ) - >rf{6) (l + %»{5) + if(2*)) . 

Applying now the Bogoliubov-Peierls and the Berezin-Lieb inequalities, we obtain: 
where 

M"M) := Tr((^-/i)P 5 )+7^(5)(l + y^ (5) + ^°(25)) . (3.22) 
Consequently, for any configuration uj one gets in the limit the estimates: 

hmsuppgV) < p(n) s$ liminf p^(/i) < limsup^M < liminf pg^ + 4 7 z/°(25)) 

+ liminf ^. (3.23) 

Z— >oo V; 

Since p^(/u) is a convex functions of /i, for any £ > we have 

Then by virtue of estimates f!3.23j) we obtain: 

limsup<9 M p"f(/i) < y(limsupp^(/x + £) - liminf p^(/z)) (3.24) 

Z — >00 * I— >oo Z— s>oo 

< i ( lip inf pj J{n + t + 47^(25)) + liminf ±-M"{l, 5, /i) - p(ji)) 

t l— >OD I— >OQ Vl 

< 7(p(/i + t + 4 7 z/°(25)) + liminf-^M a; (/,5,/i)-p( y u)) . 
Since by (I3.22p and Lemma 13.31 for any configuration u one gets: 

lim liminf —M^U, 5) = , 

54-0 Z->oo Vl 

we obtain from (I3.24p 

lim sup lim sup c^p^/i) ^ -(p(p + £) — p(aO) (3.25) 

which is valid for any t > 0. Letting £ J, leads to the second estimate (I3.18p . 

12 



Similarly one proves the first estimate (I3.17|) . Again using the convexity and (I3.23|) we 
obtain 

limsup^j^Oi) ^ i(limsupp|7(/i + t)-liminf^7(/i)) 
^ -(hmsup^(/i + 1) 



- limsupp^Gu - 4 7 z/°(25)) - liminf— AT" ^ - 4 7 ^°(26))) . 

Applying the estimate (I3.23|) to the last inequality once more, we get 

lim sup dfjffiQi) ^ I {p(ii + t)- p(ii - 4 7 z/°(25)) - lip inf ^M"(Z, 6,fji- 4rfuf(2S))) 

and in view of the limit (I3.22p and Lemma 13.31 the result follows by letting 5^0 and then, 
letting 1 1 0. □ 

Proof of Lemma 13.21 

Let 

C™* := {z E C ns : \z\ 2 < el- 
and we denote the volume of this ball by Vol(C^ s ) = ir ns ^ ns /n$r(ns)- 
We then obtain the following bound 

/ d 2 Cl ...d 2 c n ,Tr^e-^ low ^» 

= [ d 2 Cl ...d 2 c n ,Tr^e-^ low ^ Cfe » + / d 2 Cl ...d 2 c n ,Tr^e-^ low ^ Cfc » 

Jc™ s 1 JC n s\C^ s 1 

< Y^Sl maxTr ^ xe -^Mc k} ) + J_ f d 2 Cl ...d 2 c ni ( V|c fc | 2 )Tr^e-^ low ^» 



7T«* {c fc } 



Vol(C"^ 



fcG/,5 



Here we used the expectation ( ■ )i ow defined by the integrated partition function S| ou '( / u,), see 
(13. 6p . Notice that, using the explicit form of the lower symbol for the total particle number 
operator, see (13.21]) . we can represent the last integral in the form: 



d 2 Cl ...d 2 c n ,Tr^e-^ ow ^» 
1 

^maxTr^e-^ low ^» + J-tydptffdi)) f d 2 Cl . . . d 2 c„,Tr ^e~^^ c ^ 



Vol(C^ 

7T™<5 



Thus, we get that 

(1 " ^mP&W) / d 2 c, . . . d^Tr^e-^^M < (-f^) maxTr, ¥ e^W-» 
4 Jc n s y n s T(ns) / {c k } 1 
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If we put £ := aVid^™ (ji) for some a > 1, and use the Stirling formula for large rig, then 
the coefficient in the right-hand side can be estimated as: 



n*r(n*) ^ nrf- 1/2 e~^ (V^?(5)) W)+1/2 

Hence, one finally obtains the estimate: 

/ d 2 Cl ...d 2 c n4 Tr^e-^ low ^^» 

< -A^f (a^jTC/i))^ W ) V V" 1 / 2 ( I /{ ) ((y))- (,/?WV,+1/2) e^W v ' nuaVe"^^, 
1 — 1/a \ ' / {c fc } 1 

which leads to the desired result 

— ^— In f d 2 Cl . . . d 2 c n ,Tr^xe-^ low ^^» (3.26) 

< -i-lnmaxTr^e-^ low ^^» - -i-]n(l - 1/a) + ^ ln^p^)) 

+ - - ^ In (*?(*)) - — In 

□ 

Proof of Lemma 13.31 

First we note that since the projection Pg is constructed with respect to the basis of eigen- 
vectors of h®, one gets: 

v\ Vl vi 

and consequently: 

±Tr(h? - v)P S < (5-M(5)+ [ (iP l k ,v^ l k )v?(dk) 

vi J [0,6) 

= (5-M(5)+ f i$(dk)± [ dxv"(x) 

J [0,6) vi J kl 

= uf(S) ((6- (j.) + ijf dxv»(x)) . 
Thus, the ergo die theorem yields: 

limsup 1-Tr(h? - n)P s ^ v°(5)((5 - fj) + Ejv"(0) 



□ 
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4 Concluding remarks 



We discuss here the meaning of our main result and in particular the way to interpret the 
solutions of the variational problem established in Theorem 13.11 Note that the theorem 
is established for a random system in the thermodynamic limit. Its aim is to take into 
account the possibility of type III generalized condensation, which we believe is favored by 
the randomness even in interacting systems. 

First we recall a result established in [5], cf. Section 1.1. For a homogeneous system 
and the single- mode substitution at k = 0, the solution of the variational problem gives the 
total condensate density in the mode k = 0, if one adds to the Hamiltonian the zero-mode 
gauge-breaking term (quasi- average sources): 

Hi(n;rj) := Hi(p) + \[Vi (rj a + rja* ) . 

This means that after the Bogoliubov c-number substitution the solution a maXj ;(/3, p\ rj) of 
the (finite-volume) variational problem not only provides the right pressure in the thermo- 
dynamic limit, but it also coincides with quasi-averages that give the total amount of the 
condensate in the zero mode: 

lim lim a maXii (/3,/i;?7) = lim lim {a* Q a /Vi)i((3, p; rj) . 

Here (— )i((3, p; rj) is the equilibrium state defined by Hi(p;rj). 

Using a simple example, we discuss the relevance of this quasi-average approach to more 
subtle cases of the condensation of type II and III. We show that the Bogoliubov quasi- 
average sources breaking the gauge invariance [2] are able to alter the fine structure of the 
condensate reducing it to one-mode (type I). 

To see this, consider the perfect Bose-gas in a cubic three-dimensional anisotropic par- 
allelepiped A; := V^ x x V^ y x , with periodic boundary condition and a x > a y > a z , 
oix + oi y + a z = 1. The Hamiltonian is: 

fl?00 := E ( £ * ~ A*K fl fc , ( 4 -27) 

fceA* 

with quadratic spectrum and Ek=o = 0. 

It is well known that this system exhibits a generalised condensation of type II for a x = 1/2 
and of type III for a x > 1/2 for a standard critical density p c , whereas for a x < 1/2, the 
whole condensate is sitting in the mode k = , i.e, in the ground state (type I) [B]. Consider 
the system (I4.27P with the quasi-average source in a single mode k: 



#°Gu; 77) := H°(im) + y/% {f} a~ k + rj a|) . (4.28) 

Then for a fixed density p, the finite- volume equation which defines the corresponding chem- 
ical potential p = Pi(p, rj) takes the form: 

p = pi(P,p,v) := v y %2(a* k a k )?(j3, f ji,7j) (4.29) 



2 



VT ' i^e^-fO-i {e k -p) 

k^k 
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To investigate a possible condensation, we must take the thermodynamic limit in the right- 
hand side of f)4.29p . and then switch off the source, that is let \rj\ — > 0. Let us denote by 
/(/?, p) the limit of pi((3, p, r) — 0), that is the limiting density function of the gauge-invariant 
system, 

I(P,u) = lim -V af -, = / z/°(de)-5^ . 

kj^k 

with critical density p c := sup A1<0 /(/x). 

Now we have to distinguish two cases: 

(i) For any k such that lim^oo^ > 0, we obtain from (14.291) 

p = lim lim pi(/3, fi,r)) = I(/3,[i) 

r/|— s-0 I— >oo 

i.e. the quasi-average coincides with the average and we return to the analysis of the con- 
densate equation (I4.29P for r\ = 0. This gives again all possible types of condensation as a 
function of a x . 

(ii) On the other hand, if k is such that lim^oo e~ k = 0, then the condensate equation (14.291) 
yields for the quasi-average of the total particle density: 

\v\ 2 

p = lim lim p ; (/3,p, 77) = J(/3,p) + lim — - . (4.30) 

|r/|->-0 Z-»oo ??->-0 p z 

If p < Pc, then the asymptotic solution of (I4.30p is /^(p) = lim,^o h m z^oo7^(p, r]) < an d 
there is no condensation in any mode. 

If p > p c , then lim^o \rj\ 2 /~p 00 (P: r l) 2 = P ~ Pc- By explicit calculation one also gets that 
only fc-mode quasi-average is non-zero: 

1 f 1 1 Inl 2 1 

lim lim — (a* % a~ k ) l p,, rj) = lim lim <^ — + — ^ = p - p c , (4.31) 

i.e. for any the condensation is type I. Recall that the only condition on k is that the 
corresponding eigenvalue e k vanishes in the infinite volume limit. Since 

lim lim — (aSao>?(/3,A.»7) = lim lim — 1 = , (4.32) 

?7— >-0 Z->oo Vi »?— s-OZ— >oo Vi e pl - rlWU) — 1 

by virtue of (I4.31|) and (14.321) we see that the Bogoliubov quasi-average procedure not only 
transforms the generalised condensates of type II or III into a one-mode condensate (i.e., 
type I), but this mode does not even need to be the ground-state. Therefore, using the 
quasi-average approach [2J, one can force the condensate to be in any given mode k, as long 
as its energy e k vanishes in the limit I — > oo. 



Appendix 

A The approximating Hamiltonians 

In this section, we provide an explicit form of the upper and lower symbols for the Hamil- 
tonian (12.81) . For a short-hand we put Ef' 1 = Ei. We note that, as the coherent vector 
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(13. 5 p is defined as a tensor product of one-mode coherent states, its effect on each cre- 
ation/annihilation operator a\ is independent of all the others modes operators. First, we 
give an explicit form of the lower symbol of the full Hamiltonian (12. 8ft : 



+ 



+ 



kel s i>l 



+ Y CY(<Pi,4'k)(4>ki,<fri)(Ei- H))cka k i 
keh,k'ei% i^i 

+ E (53(^i»^fc)(V'fc'»0i)(- E! t - /f))4 a fc' 
k,k'ei%,k^k' i>i 



+ ^E E E <x» 



k'-qei$ 



+ ^E E E % 

fc'e/ 5 q-k+qeh 
k'-q£l% 

+ ^E E E % 

fcg/«5 fc'&Ts q:k+qel% 
k'-qeh 

+ 4EE E <«« 



fce/i fc'eia q:k+qeh 
k'-qeh 

^E E E % 

£:€/«S fc'e-ff q:k+qel% 

k'-qeii 

+ ^E E E «(« 

£:€/«S fc'S/? q.k+qeh 

k'-qei c g 

+ ^E E E % 

fcG/«S fc'6/= q:k+qel c s 
k'-qeh 

+ ^E E E % 

kels k'el% q:k+qel s 
k'-qeh 

+ wX E E % 



fcG/| fc'G/,5 q:k+qel c s 
k'-qeli 



^E E E % 

fcG/| fc'G/a q:k+qeh 

k'-geil 



CkCk 1 Clk+q a k>-q 



Ck+qCkCk' a k'-q 



Ck'-qCkCk' a 



k+q 



c k+q c k' -q c k c k' 



Ck' a k+q a k ,_ q a k 



Ck+qCk 0-k'-q a k' 



c k'-qCk Clk+q a k' 



c k+q c k'-q c k a k' 



Ck' a k +q a k'-q a k 



Ck+qCk' 0'k'~q a k 



(A.l 
(A.2 
(A.3 
(A.4 
(A.5 

(A.6 
(A.7 
(A.8 
(A.9 
(A.10 
(A.ll 
(A.12 
(A.13 
(A.14 

(A.15 



17 



+ wY E Y u(q) Ck'- q Ck' a* k+q a k (A.16) 
' keif k'el s q-.k+q&q 

k'-qelg 

+ ^E Y Y ^ 

c k+q c k'-q c k' a k 

(A.17) 

1 k£ll k'els q:k+q£l s 
k'-qels 

+ ^Y Y Y ""fa) a* k+q a* k ,_ q a k ,a k (A.18) 

' fee/| k'&II q:k+q&I c s 
k'-qel% 



+ ^yYY Y u(q) c k + q a* k ,_ q a k ,a k (A.19) 

1 kel$ k'£ll q:k+q£lg 

k'-qeq 

+ ^jE E Y u(q)ck>- q a* k+q a k >ak (A.20) 

1 fcei| fc'eij q-k+q£l c & 
k'-qels 

+ ^yYY Y u(q)c k+q c kl - q a kl a k (A.21) 



' fceij fc'eij q-.k+qeig 
k'-qeis 

Now, we give an explicit form of the upper symbols. We recall the general form of this 
symbols for polynomials in the creation/annihilation operators of some mode k G Is'- 

(a k ) u P = c k , (al) up = c k , (a k a k ) u P = (c k ) 2 , (a* k a k ) up = (c k ) 2 
(a* k a k ) up = \c k \ 2 - 1, (a* k a* k a k a k ) up = |c fc | 4 - 4|c fc | 2 + 2 

Note that, since the interaction term of the Hamiltonian term considered on its own does have 
a momentum conservation law, it is not possible to get exactly three out of four operators 
in the same mode k. In view of this, it can be seen that the lower and upper symbols 
of the Hamiltonian will differ only when two or four operators in the same mode appears, 
that is only terms (TO) . (TA~9|) . (EQ2]) . flA~T3|) . flAl5j) . flAl6j) and (1X17)1 differs in both 
approximating Hamiltonians. 

Splitting further the sums in these terms leads finally to the final upper symbol of the 
Hamiltonian 

H^,{c k }) = H l r(fi,{c k }) + ^,{c k }), 

where 

k(Mc*}) = J2{^2\(Mk)\ 2 (Ei-^) (A.22) 

+ -^( 2 E^°) + E ^(°) + E E 

' fee/ 5 fe,fc'G/«5 g^o kei s m-q 

k^k' 

+ ^7( 4 E^°)I^I 2 + E ^(o)(M 2 + M 2 ) + ]T £ %)(| Cfc | 



Cfc' 



2- 



fcg/,5 k,k'el s q¥=0 k€l s nl-q 

k^k' k'ei s m-q 

k'~k=q 

77 ( 2 E "(°) E a fc' afc ' + E E u(q)a* k a k + «(?K< 



2V, 



i k a k j . 

keig fc'eij g^o fce/ 5 g^o fceijni-q 

fc'eijni+q k'eh 

k'—k=q k'—k=q 
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